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No magnetic field variation with pulsation phase 
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Abstract. We present an analysis of 210 high-resolution time-resolved spectropolarimetric observations of the roAp star yEqu 
obtained over three nights in August and September 2003. Radial velocity variations due to p-mode non-radial pulsations are 
clearly detected in the lines of rare-earth elements, in particular Pr in, Nd n and Nd in. In contrast, we find absolutely no evidence 
for the variation of the mean longitudinal magnetic field over the pulsation period in y Equ at the level of 1 10-240 G which was 
recently reported by Leone & Kurtz 1 2003 1. Our investigation of the variability of circularly polarized profiles of 13 Nd in lines 
demonstrates that, at the 3cr confidence level, no magnetic field variation with an amplitude above a 40-60 G was present in 
y Equ during our monitoring of this star. 
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1. Introduction 

Since the discovery of high-overtone /?-mode pulsations in the 
rapidly oscillating Ap (roAp) stars it became clear that, un- 
like in other stellar pulsators, strong magnetic fields in these 
stars have a defining role in exciting the oscillations and shap- 
ing main pulsational properties. It was found that the ampli- 
tude and phase of the rapid light variation are modulated by 
the stellar rotation and that the phases of the extrema of mag- 
netic field and pulsational amplitude typically coincide with 
each other. A phenomenological oblique pulsator model (Kurtz 
1982 1 attributed main characteristics of the photometric pul- 
sational variation of roAp stars to the oblique dipole modes, 
aligned with the axis of dipolar magnetic field. 

The apparent connection between global magnetic field of 
roAp stars and their non-radial pulsations led Hubrig et al. 
( 2003 1 to propose that atmospheric magnetic structure may also 
change over the pulsation period and hence pulsational vari- 
ability of the longitudinal magnetic field can eventually be ob- 
served. Hubrig et al. l|5p03 1 searched for the magnetic varia- 
tions in a sample of six roAp stars but failed to detect any co- 
herent magnetic variability in low-resolution Zeeman spectra. 
This null result is not unexpected given that in roAp stars lines 
of different elements display dramatically different amplitudes 
and phases of oscillations, which are combined destructively 
and result in an extremely weak pulsational signal when aver- 
aged over broad wavelength regions. 
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High-resolution analysis of Stokes profiles of individual 
metal lines is a much more promising tool for mapping atmo- 
spheric pulsational disturbances and investigating possible os- 
cillations of magnetic field. In a recent paper Leone & Kurtz 
(2003 hereafter LK) announced the discovery of variation of 
the mean longitudinal field (B z ) with the 12. 1-minute pulsation 
period in one of the brightest roAp star y Equ. Magnetic vari- 
ability with an amplitude between 110 and 240 G was inferred 
from the study of 4 Ndra lines in only 18 circularly polarized 
spectra of yEqu, which makes the alleged detection of mag- 
netic variation a statistically marginal result. With the aim of 
verifying the existence of rapid magnetic variability in yEqu 
we decided to acquire and analyse a much more extensive time 
series spectropolarimetric observations of this star. 

2. Observations 

We obtained time-resolved circularly polarized spectra of 
y Equ using the cross-dispersed Nasmyth Echelle Spectrometer 
(NES, Panchuk et al. [1999 1 installed at the 6-m telescope of 
the Russian Special Astrophysical Observatory. The NES in- 
strument is equipped with a Zeeman analyzer consisting of a 
quarter-wave and half-wave retarders and a calcite plate sepa- 
rating the beams with opposite circular polarization. A 2Kx2K 
Loral CCD detector allowed us to record 26 echelle orders with 
nearly complete wavelength coverage in the 4520-6000 A re- 
gion with a spectral resolving power of R ss 38 000. 

On each of the nights of August 19 and 20, 2003 we ob- 
served yEqu for 2.4 h and have obtained continuous sequences 



2 



O. Kochukhov et al.: No magnetic field variation with pulsation phase in the roAp star y Equulei 



of 70 Stokes / and V spectra using 80 s integration time with a 
readout time of 42 s . An additional series of 70 polarized spectra 
of yEqu was secured on September 11, 2003, when exposure 
times of 120 s were used. On each of the observing nights the 
peak S IN ~ 60-80 was reached in individual time-resolved 
spectra. The software package REDUCE (Piskunov & Valenti 
2002 1 was used for the optimal extraction of the NES spectra. A 
wavelength scale with an internal accuracy of 50-70 m s was 
established by means of a 2-D wavelength calibration proce- 
dure, which made use of ^700 ThAr lines in all echelle orders. 
The left- and right-hand circularly polarized (LCP and RCP) 
beams were calibrated independently using the corresponding 
beams of the reference spectrum taken before and after stellar 
observations. 

3. Identification of Ndm lines 

Investigation of pulsational variability in the spectral lines of 
the rare-earth elements (REE) offers unique potential to ob- 
serve non-radial /9-mode oscillations and to study in detail 
the structure of pulsating cavity in roAp stars. It was recently 
demonstrated by Kochukhov & Ryabchikova (2001a, 2001b) 
that the largest pulsational variability (RV amplitudes up to 
1 kms~') in individual metal lines in yEqu and other roAp 
stars is invariably found in the REE spectral lines, in partic- 
ular strong absorption features of doubly ionized Pr and Nd, 
while pulsations could not be detected in the lines of light and 
iron-peak elements. Since the lines of each REE ion show very 
similar pulsational behaviour, the accuracy of our search for 
magnetic oscillations in yEqu can be substantially improved 
by simultaneous multiline analysis of Nd m lines, which consti- 
tute the most numerous group of strong lines with pronounced 
variability in the spectrum region available in our observations. 

In the present study magnetic field and RV variations 
were measured using 13 unblended Ndm lines some of which 
were previously unclassified, but appeared in laboratory spec- 
tra (Crosswhite ll976l Aldenius 2001 1 as well as in the spectra 
of roAp stars known to have Nd m lines of abnormal intensity 
(Ryabchikova et al. l2001> . We emphasize that all newly classi- 
fied Nd in lines show strong pulsational signatures, consistent 
with the behaviour of known lines of Nd m, which by itself pro- 
vides a convincing support for our line identification. Table [2 
contains wavelengths of Nd m lines, proposed transition iden- 
tification based on the energy level calculations by Zhang et 
al. ([2002 1 and the effective Lande factors which were calcu- 
lated from the ^-factors published by Bord (2000) and Zhang 
et al. (2002). The details of the line classification will be pre- 
sented elsewhere, here we note that the adopted Lande factors 
and Zeeman patterns are consistent with the observed Zeeman 
splitting in the sharp-lined Ap star HD 144897 with the surface 
magnetic field of w 9 kG. 

4. Radial velocity and magnetic field 
measurements 

The pulsational analysis of y Equ followed the procedure out- 
lined in Kochukhov & Ryabchikova (2001a). For each spec- 
tral line in the available wavelength region we determined the 



Table 1. Nd in lines used in the RV and magnetic field analysis 
of yEqu. The columns give laboratory wavelength, effective 
Lande factor, excitation potentials and J values of the lower 
and upper levels. 
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"Previously identified and classified Nd m lines 

"Lines included in Crosswhite's 1 1976) unpublished list 

'Lines detected in the laboratory spectrum (Aldenius 2001 and private 

communication) 

"Reclassified and newly classified lines 

center-of-gravity as a function of time in the Stokes / spectra. 
Then a refined analysis of polarization spectra was carried out 
for the group of Nd m lines aimed at detection of the magnetic 
variability. An instantaneous position of the line center was de- 
termined consistently for the RCP and LCP profiles by integrat- 
ing over the spectral points below certain residual depth. The 
same vertical cutoff was applied to the RCP and LCP spectra, 
but it had to be chosen individually for each line due to differ- 
ent blending in the far wings. This results in up to a factor of 
two difference in the RV amplitudes for studied Nd m lines, al- 
though depth effects may also be partially responsible for that. 
The RV and (B z ) were derived from the center-of-gravities of 
the RCP and LCP profiles, Ar and A^, using the expressions 

(B ) = + A L 

2.34x10- 1 W<^r> + <^l» 2 ' m 

where c is the speed of light and (/1r,l) correspond to the time- 
averaged line center positions. 

Preliminary time series analysis of the Nd m lines showed 
negligible difference between the periods and phases of their 
RV variability. Therefore, in subsequent least squares sinu- 
soidal fits we allowed for individual amplitudes, but used the 
same phase and pulsation period for all Nd m lines, effectively 
coadding information from many lines into a single high S /N 
RV measure. The same strategy was used to construct the av- 
erage variation of the longitudinal field after subtracting a con- 
stant offset from (B z ) measurements obtained for each of the 
Ndm lines. Results of our investigation of the neodimium line 
variability are summarized in Table ^ an d are illustrated in 
Figs. H an d 13 which show the average RV curves (scaled to 
nightly mean amplitudes) as well as the average variation of 
(B z ) and corresponding amplitude spectra. 
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Fig. 1. Pulsational variation of the radial velocity and longitudinal magnetic field measured for Ndm lines on the nights of 
August 19 and 20, 2003. The average pulsational RV curves a) are folded with the best-fit oscillation period and are compared 
with the average variation of (B z ) c). Panels b) and d) illustrate the amplitude spectra for the RV and longitudinal magnetic field. 
The vertical dashed lines show photometric pulsational periods of y Equ. The horizontal dashed line in panels d) indicates the 
maximum amplitude of pulsational {B z ) variation reported by LK. 



Table 2. Results of our analysis of the time-resolved polarized spectra of yEqu. The columns list UT date of observation, 
reference Julian date, number of analysed spectra and a least-squares estimate of the pulsational period, average amplitude and 
phase (measured in fractions of the period) of the pulsational RV variation in 13 Ndm lines. The last three columns report the 
amplitude of the magnetic variation, standard deviation of the (B z ) measurements and an amplitude which would have been 
detected with out data at the 3cr confidence level. 
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For all three nights the pulsation period detected in our data 
is consistent with one of the frequencies (P = 12.201 min) 
identified in photometry (Martinez et al. 1996 1. Despite unam- 
biguous detection of the pulsational RV shifts, no evidence for 
magnetic variability is found. A rigorous statistical analysis in- 
dicates that, at the 3<x confidence level, no magnetic variability 
with A{B Z ) > 40-60 G is seen in y Equ during the nights of our 
observations. 

In addition to the study of the Nd m lines, we identified a 
number of Prm and Ndn features showing clear pulsational 
variability. These lines are not so numerous and strong as to 
provide an independent and accurate diagnostic of possible 
rapid magnetic changes. Nevertheless, analysis of these ions is 
consistent with the results obtained using the Ndm lines: two 
lines of Prm with the average RV amplitude between 950 m s _1 
(19 Aug) and 280 ms" 1 (11 Sep) show A(B Z ) < 100 G, while 
A(B Z ) w 60+50 G is derived from four Nd 11 lines. Remarkably, 
the lines of singly ionized Nd, which could not be investi- 
gated in our first spectroscopic study of pulsations in yEqu 
(Kochukhov & Ryabchikova 200 1 a 1, show very high average 
RV amplitudes of 1416 ms" 1 (19 Aug), 1118 ms" 1 (20 Aug) 



and414ms~ 1 (11 Sep) and a constant phase shift of 0.13+0.02 
of pulsation period with respect to the RV variation of Nd m. 

The procedure of the magnetic and RV measurements ap- 
plied to REE lines was also tested on a sample of 19 unblended 
strong lines of Ca, Cr and Fe, which are not expected to show 
any pulsational variability of either RV or (B z ). From the Aug 
19 data we obtained RV amplitude of 30 + 19 ms~', A(B Z ) = 
18+14Gand standard deviation cr((B z )) = 84Gfor these lines, 
in good quantitative agreement with the REE results. This pro- 
vides a complementory indirect argument for the absence of 
magnetic variability in REE lines. Since the variation of mag- 
netic field is predicted to be proportional to the RV oscillations 
(see Hubrig et al. 2003 1, even for marginal magnetic variation 
one would expect to find ct((B z )ree) » cr((B z )c a ,cr,Fe), but this 
is clearly not observed in yEqu. 

5. Discussion 

What may be the reason for a major discrepancy between the 
results reported by LK and the outcome of our study of yEqu? 
We consider it improbable that the star itself has changed 
its pulsational behaviour dramatically. Theoretical estimate of 



O. Kochukhov et al.: No magnetic field variation with pulsation phase in the roAp star y Equulei 



11/09/2003 



Nd III 5802.54 A 



Nd III 5845.07 A 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Pulsation phase 



1 200 



v -200 
< 

-400 : 
-600 ': 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Pulsation phase 




Fig. 2. The same as Fig.[Qfor the spectropolarimetric observa- 
tions of y Equ obtained on the night of September 11, 2003. 



possible magnetic variability (Hubrig et al. 2003 1 hints that 
magnetic and velocity amplitudes are proportional to each 
other. Therefore, if magnetic variability really exists, on the 
nights of Aug 19 and 20 (when RV amplitudes of Ndm lines 
were about twice the amplitude measured by LK) we should 
have observed even larger (B z ) variation than found by LK. 

It is possible that the difference between our null result 
and the alleged positive detection of the magnetic variability 
by LK is related to different approaches to measure (B z ) and 
the effects of blending of the Nd in lines studied by LK. We 
derive longitudinal field from consistent center-of-gravity de- 
termination in the RCP and LCP profiles, while LK analysed 
the first moment of the Stokes V spectra. Strong Ndm lines 
tend to show extended wings and have unusually wide Stokes 
V signatures. Hence, an insignificant blending in Stokes / may 
correspond to a substantial overlap of the Stokes V profiles of 
neighbouring lines. The Nd m 5845.07 A, which LK offered as 
the most prominent case of (B z ) variation, represents a clear ex- 
ample of the blending problem. This line is blended by the vari- 
able line of Prm at A 5844.41 A, which appears to be reason- 
ably well separated from the Nd m line in Stokes /, but signifi- 
cantly distorts the corresponding Stokes V signature. In Fig. [3] 
we use our Aug 19 data to illustrate that the blending of the 
Ndm 5845.07 A results in spurious variation of (B z ) when in- 
tegrating Stokes V signal for this line within a fixed wavelength 
range. In contrast, no variation is seen in the nearby line of 
Ndm 5802.54 A which is free from any blends and has higher 
magnetic sensitivity. 

We also point out that the results of the RV and magnetic 
measurements presented by LK for individual Nd m lines are 
inconsistent. The reported phase shift between (B z ) curves of, 
e.g., the Ndm 5845.07 and 6145.07 A, reaches 0.20 of the pul- 
sation period (with a typically phase error not exceeding 0.03), 
while no such shift is seen in RV variation. Unfortunately, LK 
do not comment on this peculiar behaviour of different Nd m 
lines. 
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Fig. 3. The average Stokes / a) and V b) profiles of the Nd m 
lines at A 5802.54 and 5845.07 A. The vertical dashed lines 
illustrate selection of the spectral region for the RV and mag- 
netic measurements. The panel c) shows variation of the (B z ) 
measure obtained from the Stokes V profiles. Blending by the 
Prm 5844.41 A leads to spurious variation of {B z ) determined 
for the Ndm 5845.07 A. 

Based on our observations we conclude that we see no evi- 
dence for pulsational variation of magnetic field in y Equ con- 
trary to the recent claim by LK. We show that systematic prob- 
lems exist in the latter study and hence the alleged detection of 
the magnetic variability in y Equ may be spurious. 
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